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O^' Abstract 

Generic channel and ridge structures are known to appear in the 
KN ' time-dependent position probabiUty distribution of a one-dimensional 

j^ ■ quantum particle confined to a box. These structures are shown to 

have a detailed quantitative explanation in terms of a travelling-wave 
decomposition of the probability density, wherein each contributing 
term corresponds simultaneously to (i) a real wave propagating at 
a quantised velocity and (ii) to the time-averaged structure of the 
position distribution along a quantised direction in spacetime. The 
approach leads to new predictions of channel locations, widths and 
depths, and is able to provide more structural details than earlier ap- 
proaches based on partial interference and Wigner functions. Results 
are also applicable to light diffracted by a periodic grating, and to the 
quantum rigid rotator. 



1 Introduction 

The position probability distribution P{x, t) of a one- dimensional quantum 
particle may be represented as a probability landscape in spacetime, where 
hills and valleys correspond to regions of high and low probability density 
respectively. It has recently been discovered that generic two-dimensional 
structures appear across this probability landscape for a wide variety of po- 
tentials and initial conditions [|l|]-|0. The highly-patterned nature of these 
structures has led to them being called "quantum carpets" . 

For the simplest case, of a one-dimensional particle moving freely between 
two end walls, the carpet patterns are linear, with quantised slopes and inter- 
cepts [0 ^, 3 More generally, however, the patterns are curved for particles 
moving under the influence of a potential @, H, P. Colour plots of quantum 
carpets, for a range of examples, may be found in [Q and 0. 

Analogous carpet structures also arise for light diffracted by a one-dimensional 
periodic grating in the paraxial approximation 0, §]. In this case the time 
variable t is replaced by the propagation distance z of the diffracted beam, 
and P{x, z) is the intensity distribution of the beam. Further, it will be 
seen that a generalised form of these carpet structures arises for the angular 
distribution P(0, t) of the two-dimensional quantum rigid rotator. 

Approaches based on interference between energy amplitudes p, |^ , Wigner 
functions @], and Greens function degeneracies P, |] have been used to ex- 
plain the observed carpet patterns for various cases, and in particular the 
observed quantisation of slopes of the linear structures for the confined one- 
dimensional particle is well understood. However, even in this relatively 
simple case a number of generic features have not yet been provided with a 
general explanation, including the depths and widths of channel structures, 
the decrease of the latter two quantities with slope, and the "chopped" nature 
of ridge structures (see Fig. 1). 

It will be shown in section 2 that such features can be understood to 
a limited extent by generalising an approach first developed by Berry [Q], 
based on the destructive interference of energy amplitudes (see also 0]). 
This generalised destructive interference approach in particular predicts val- 
ley structures along the lines 

X = kVt + IL, kl = even integer (1) 

in spacetime for a particle confined to a one-dimensional box with end walls at 



X = and x = L, providing that the energy amphtude differences lipn — i'n+kl 
of the wavefunction are sufficiently small for all n. Here k and / are integers, 
and 

V = Trh/{2ML) (2) 

is a natural speed defined via the particle mass M. 

The destructive interference approach can be further used to show that 
the average channel depth is expected to decrease as k increases (important in 
resolving the Olbers-type paradox that any channels at all can be discerned, 
given that an infinite number of them are predicted by (|l]) in any given space- 
time region). Moreover, as shown in the Appendix, it may be generalised to 
predict the locations of channel and ridge structures for particles moving in 
arbitrary one-dimensional potentials, with results in agreement with [Q . 

However, the above approach is not able to explain a number of other 
observed features. First, it is completely silent on the location of structures 
when the condition of small energy amplitude differences is not met. Second, 
as will be shown by example, it fails to predict all observed channel and ridge 
structures even when this condition is met. Finally, the approach yields no 
information on the shapes and widths of the various structures. A new 
approach is thus clearly called for. 

Such an approach is developed in section 3, based on a representation 
of the probability landscape as a superposition of travelling waves with ve- 
locities which are integer multiples of V, propagating against a constant 
background. The primary usefulness of this decomposition is that the aver- 
age of the probability distribution along any given direction in spacetime (eg, 
corresponding to a channel or ridge) is described by at most a single one of 
these travelling waves. In particular, the wave corresponding to velocity kV 
provides specific, exact information on the average locations, shapes, depths 
and widths of all channel and ridge structures with slope kV. Connections 
with the destructive interference approach are briefly discussed. 

The averaging of P{x,t) along particular spacetime directions, to pre- 
dict locations of corresponding channel structures, was first (independently) 
suggested by Berry and Bodenschatz p in the optical grating context, who 
were also motivated by limitations of the destructive interference approach. 
That these averages arise naturally from a travelling wave decomposition of 
P{x, t), and can hence be simply superposed to recover the complete probabil- 
ity landscape (and to approximate it with arbitrary accuracy by considering 



only waves up to a given maximum speed) has not previously been realised. 

In section 4 it is shown further that each travelling wave can be written 
as the sum of two Wigner functions plus an energy amplitude term. This 
form is useful for showing how localisation and symmetry properties of the 
initial probability distribution can enhance or suppress carpet patterns. In 
Sec. 5 predictions are verified via the examples of a uniform initial wavefunc- 
tion and an approximately Gaussian initial wavefunction, for which analytic 
expressions are obtained for the travelling waves. 

Generalisations of the results to periodic optical gratings and to the quan- 
tum rigid rotator are given in section 6, with conclusions presented in section 
7. 

2 Destructive interference approach 

2.1 Derivation 

Berry derived the locations of channel structures of a one- dimensional particle 
moving freely between two end walls, for the case of a uniform initial wave- 
function [^ . Here his derivation is generalised to predict the quantum carpet 
structure for a broad class of initial wavef unctions, and in the Appendix it 
is shown how this approach may be further generalised to predict the curved 
carpet structures corresponding to arbitrary one-dimensional potentials. 

Schrodinger's equation for a free particle of mass M confined to the in- 
terval (0,L), under the usual boundary conditions that the wavefunction 
vanishes at the endpoints, has the general solution 

oo 

V'(x,t) = (2/L)i/2^7/;„e-'""'^*/^sinn7rx/L (3) 

n=l 

r oo oo 

= -1/V2L ^7/'„e^'^+(^'*'")-^V'„e^'^-(^'*'") , (4) 

.n=l 71=1 

where ipn denotes the n-th energy amphtude, V is the speed defined in (0), 
and the phases 0± are given by 

0±(x, t,n) = (±nx — n^V"t)7r/L. (5) 

The position amplitude ip{x,t), and hence the position probability dis- 
tribution P{x,t) = |-?/'(a;, t)p, will be small if partial cancellation can be 



arranged between the two summations in (^. This is possible in particular 
if the term corresponding to summation index n in the first summation par- 
tially cancels with the term corresponding to summation index n + k m the 
second summation, where k is some fixed integer. Indeed, it is not difficult 
to show that (1) can be rearranged (to within an overall phase factor) as 

oo 

n=l 

\k\ 

+ (2L)-i/2^7/;„e'^±("'*''^), (6) 

n=l 

where the upper (lower) phase subscript is chosen when k is positive (nega- 
tive). Thus significant cancellation between energy amplitudes can take place 
in the first summation, leading to a small value of P(x, t), providing that the 
phase-matching condition 

0+(a:, t,n -|- fc) = 0_(x, t,n) (mod 27r) (7) 

is met for all n. 

It is the phase condition (|^) which leads to (|1|). In particular, note that 
channel structures in the probability landscape correspond to spacetime tra- 
jectories x{t) along which the position probability density is relatively low, 
and hence in particular to trajectories for which the phase condition (0) is 
satisfied at all points. Substituting (|) into (|^, differentiating with respect to 
time, and comparing consecutive values of ra, yields the condition x{t) = kV 
for the slope of such a trajectory. Substituting t = in (^ further yields, on 
comparison for consecutive values of n, condition (|1|). 

Note that while condition (||) thus follows from destructive interference 
between pairs of energy amplitudes in (|^), two further conditions are also 
necessary in general for channel structures to be observed. In particular, 
substitution of ([1|) in (^ gives the expression 

oo 

i,{kVt + mt) = (2L)-V2^(_l)-'(^„^|,|-^„)e— ('^+"=1)-^*/^ 

n=l 

\k\ 

+ (2L)-^/2^(-l)"Vne-*"^"-l'l)"^*/^ (8) 

n=l 



for the position amplitude along the trajectory. Thus for the trajectory to 
correspond to a channel one requires further that (i) the differences {ipn — 
ipn+\k\\ Ojre small , and (ii) \k\ itself is sufficiently small. These requirements 
guarantee that the first and second summations in (H) yield, respectively, 
relatively small contributions to the total amplitude, and hence to P{x,t). 
They hold for a wide group of initial wavef unctions, as discussed in section 
2.3, and hence condition (P has a wide predictive power. 

2.2 Channel depths and ridge heights 

Equation (H) leads to a simple estimate of channel depth. In particular, 
the quadratic n-dependence of the phases in (Q) implies that they are quasi- 
random in time, and hence the average probabilty density along a given 
channel can be estimated as 

oo \k\ 

|2 



P{kVt + lL,t) ^ {2L)~'J2\^n-^Pn+\k\\^+{'^L)-'J2\^r. 

n=l n=l 

oo 

= L-i(l-^Re{CV^„+|fc|}). (9) 



n=l 



Note that for slowly varying energy amplitudes this expression will typically 
increase as |fc| increases, i.e., deep channels correspond to small values of \k\. 
It will be seen in section 3 that this expression is close to the exact average 
depth of the channel. 

A particularly simple example is when the initial wavefunction is an 
equally weighted superposition of N consecutive energy eigenstates, i.e., 
ijjn = 1/VN ioT M < n < M + N and ijjn = otherwise, for some M > 0. 
From (Q) it follows that relatively good destructive interference takes place 
for I A; I '^ N, and from (^ that the average density along a channel is approx- 
imately given by |A;|/(A^L). Thus deep channels correspond to small values 
of \k\, as expected. Note moreover from (§) that destructive interference is 
impossible for \k\ > N, and hence only a finite number of channels can be 
observed in this case. A previous analytic study of this example for M = 
confirms these results [^. 

Constructive interference between energy amplitudes in the first sum- 
mation in (^ corresponds to adding tt to one side of the phase-matching 



condition (|^). It follows that constructive interference takes place along tra- 
jectories as per (|l]), except that the product kl is now restricted to be odd. 
Thus "ridges" are predicted along such lines in the probability landscape, in 
agreement with numerical observations [0] - [Q, 0. Moreover, the average 
height of these ridges may be estimated similarly to channel depths above, 
yielding an expression similar to (^ but with subtraction replaced by addi- 
tion. Thus for the above example the average ridge height is predicted to be 
{2N-\k\)/{NL) for \k\ < iV. 

Finally, note that fluctuations of channel depth along a given channel 
are constrained to be small, simply because they are bounded by the rela- 
tively small average probability density (which must remain positive along 
the channel). In contrast, fluctuations of ridge height are only bounded by 
the relatively large average ridge height, and indeed are expected to be of 
the same order from the quasi-random nature of the phases appearing in the 
analogue of (^ for constructive interference (where the minus sign in the first 
summation is replaced by a plus sign). It is these relatively large fluctuations 
which lead to the observed generic "chopped" nature of ridge structures . 

2.3 Examples 

The usefulness and limitations of the generalised destructive interference ap- 
proach is investigated here via two generic examples. It will be shown in 
particular that channels are expected as per ([^) for well-localised initial 
wavepackets; and conversely that destructive interference fails to account 
for all observed channels in the case of wavepackets with periodically spaced 
energy amplitudes. 

Example (i) Momentum amplitudes and localised wavepackets: 
The probability landscapes of wavefunctions initially well-localised in a one- 
dimensional box have been studied in a number of special cases [^, ^, ^ 
p. Here it is shown that the destructive interference approach provides a 
general explanation of channel locations for all initial wavepackets which have 
a slowly varying momentum amplitude distribution, and for well-localised 
wavepackets in particular. 

Let ipip) denote the momentum amplitude distribution of the initial 
wavepacket, i.e., 

i){p) = {2^h)-^'^ e-'P"/V(a;, 0)rfa:. (10) 

J —oo 



It follows from (||) and (]I0|), recalling that ipi^x, 0) vanishes outside the inter- 
val (0, L), that the energy amplitudes ipn can be expressed as 



L 



ipn = (2/L)^/2 / sm{mix/L)ij{x,0)dx (11) 





= i{7rh/Ly/^[ip{mTn/L) --^{-nTth/L)]. (12) 

From (|^) and (p!^ ) it follows the destructive interference approach predicts 
channels of slope kV, as per (||), providing that the momentum amplitude 
distribution ip{p) varies sufficiently slowly over any range of length \k\TTh/L. 
In particular, if the initial wavepacket V^(x, 0) is smooth and well-localised 
within the box then the momentum amplitude distribution will typically be 
broad, with root mean square variance Ap say. Thus channels of slope kV 
are predicted for 

\k\Trh/L < Ap. (13) 

Noting that the Heisenberg uncertainty relation implies that Ap > h/{2Ax), 
one in particular expects to see channels of slope kV for all smooth localised 
wavepackets such that |/!;| -C L/{2'kAx). These general predictions are well 
borne out by the approximately Gaussian initial wavepackets studied previ- 
ouslyi, i-i. 

Example (ii) Periodically spaced energy amplitudes: The second 
example to be considered here is the case where the non-zero energy ampli- 
tudes of the wavefunction are periodically spaced, i.e., 

'?/'„ = for n ^ r mod p (14) 

for some period p and fixed integer r. 

The case p = 1 is trivial, with channels predicted as per (|I|) just as 
before. The case p = 2 corresponds to initial wavefunctions which are either 
symmetric or antisymmetric about x = L/2, as r is odd or even respectively, 
as can be directly seen from the form of the wavefunction in (|^). The case 
p = 3 will be seen to be of special significance, providing an example where 
the destructive interference approach breaks down. 

From ( |T4|) it is seen that amplitudes ipn and ipn+\k\ can only destructively 
interfere in the first summation in (^ when A; is a multiple of the period p 
and n = r mod p. Accordingly, substituting pk for k and pn + r for n in the 
phase-matching condition (|^ leads to the modification 

X = kpVt + IL/p, {k + 2r/p)l = even integer (15) 



of (|1|), for trajectories which correspond to channels in the probabihty land- 
scape. Condition (|l^) is equivalent to (|l]) for p = 1. More generally the 
predicted channel locations depend upon both p and r, have slopes which 
are multiples oi pV, and intersect the x-axis at multiples of L/p. 

For example, for p = 2 and r = 1 (symmetric initial wavefunctions) , it 
follows that channels correspond to the trajectories x = 2kVt + lL/2, such 
that {k + l)l is even. This prediction is equivalent to Eq. (52) of §, obtained 
there for the special case of an initially uniform wavefunction. 

The case p = 3 provides an example which demonstrates an incomplete- 
ness of the destructive interference approach. In particular, consider an initial 
wavefunction given by the equally-weighted superposition 

N-l 

^(x, 0) = Ar-^/2 J2 sin[(3n + l)nx/L]. (16) 

n=0 

This corresponds to p = 3 and r = 1 in (|T^), and hence from (|T^ no 
channels starting from x = L/3 are predicted. However, the density plot of 
the corresponding probability landscape for this example, shown in Figure 1 
for A^ = 20, shows that a channel of slope —3V is in fact associated with this 
starting point. The appearance of this unexpected channel will be explained 
in the following section. 

3 Travelling wave approach 

3.1 Structure function decomposition 

It has been seen that the destructive interference approach successfully pre- 
dicts a number of generic properties of quantum carpets, including the loca- 
tions of channel and ridge structures and their corresponding average depths 
and heights. It does not, however, yield information on the shapes and widths 
of these structures; does not predict a// observed structures (Figure 1); and is 
in any case limited in applicability to wavefunctions with energy amplitudes 
which are slowly varying over ranges of length \k\. Thus a new, more general 
approach is desirable. 

To introduce such an approach, it is first convenient to rewrite the wave- 



function in (|^) in the form |^ 

oo 

^{x,t) = -i{2L)-^/^ ^ ^,e'(""-"'^*)"/^, (17) 

n=— oo 

where one extends the energy amphtude coefficients ^^i to negative values of 
n via the definition 

1p-n = -V'n- (18) 

From ([T7|) one immediately has 

oo 

P{x,t) = (21)-^ J2 ^;;7/'„e-*("-")[^-("^+")^*]"/-^ (19) 

m,n=— oo 

for the position probability distribution. It may be seen that each term in 
this summation either has no spacetime dependence {m = n), or is a plane 
wave with velocity {m + n)V {m ^ n). Collecting the m = n terms into 
a constant background term, and grouping terms corresponding to waves 
propagating at velocity kV, one obtains the decomposition 



P{x,t) = L-^ 



- fx-kVt' 

1+ E s.[-^— 

k=—oo \ / 



(20) 



of the probability landscape, where the "structure functions" Sk are given by 

1 °° 1 

Sk{z) = ^ E rm^k-me-'^'-^-'^-'-^l^kM'. (21) 

m=— oo 

Here ipk/2 is defined to be zero when k is odd. 

Equation ( PU[ ) is the travelling wave decomposition referred to in the 
Introduction. The first term, L~^, is a constant background term for the 
probability distribution, and the structure function Sk is a real travelling 
wave travelling at velocity kV. Note that S^ is defined on the entire real 
axis, and satisfies the relations 

Skiz + 1) = i~l)''Sk{z), S^kiz) = Sk{-z), (22) 

where the latter of these follows using (p!8|). 

The travelling wave decomposition (pOD provides a physical picture for the 
generation of the probability landscape P{x, t), as the superposition of waves 

10 



of discrete velocities propagating in spacetinie against a constant background 
probability. However, the primary usefulness of this decomposition arises 
from its relationship to the statistical properties of the probability landscape. 



In particular, from (pi]) , the time average of Sk{z + at) is equal to Sk{z) if 
a = and vanishes otherwise. It follows via (^) that the average oi P{x, t) 
along the linear trajectory x{t) = Xq + kVt is given by 



P{xo + kVt,t)=L-^[l + Su{xo/L)] (23) 

for any integer k, and by L~^ when k is not an integer. Thus the average 
distribution along any direction in spacetime involves at most one structure 
function. 



3.2 Locations, shapes and widths of Unear structures 



From ( P3D it is seen that Sk{z) contains all information pertaining to the 
average properties of linear structures of slope kV in a quantum carpet. In 
particular, channel/ridge structures are associated with those trajectories 
X = xo + kVt for which xq/L corresponds to minima/maxima respectively of 
Sk- Moreover, the average shape of these stuctures corresponds to the shape 
of the associated structure function. 

For example, 1 + S^^{z) is plotted in Figure 2 for the wavefunction (|1^) 
(with A^ = 20). The deep maximum in Figure 2 in the vicinity oi z = 
1/3 implies, via relation (^), that a well-defined channel structure of slope 
—31/ crosses the x-axis at L/3. This channel is precisely the unexpected 
channel observed in Figure 1, which was not predicted by the destructive 
interference approach of section 2. Note also from Figure 2 that a ridge 
structure immediately parallel to the left of this channel is also predicted, as 
indeed may also be observed in Figure 1. Since the structure functions in (pl|) 
can be trivially evaluated as sums of geometric series for this example, the 
average properties of all linear structures can in fact be calculated analytically 
via ( P5D if desired. 

In Figure 3 the travelling wave decomposition (EH) is directly illustrated, 



again for the wavefunction ([16D with A^ = 20, where only terms with |fc| < 5 
have been included in the summation. Comparison with Figure 1 demon- 
strates that all linear carpet structures with slope less than or equal to bV 
in magnitude are reproduced in Figure 3. 
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Note finally that one can make use of (^3]) to define the average horizontal 

width of a given channel or ridge structure. In particular, suppose such a 

structure corresponds to a trajectory passing through the x-axis at Xq. Define 

x+ and x_ to be the points to the right and left respectively of Xq for which 

the average probability distribution first becomes equal to the background 

probability L~^. The difference x+ — a;_ is then a natural measure of the 

horizontal width of the structure. From ( p3|) these points correspond to the 

zeroes of the structure function lying either side of xq/L, where xq/L itself 

corresponds to a minimum or maximum of the structure function. Note for a 

structure of slope kV that a corresponding measure of width in the direction 

perpendicular to the structure is obtained by dividing the horizontal width 
by (l + fcV2)V2. 

3.3 Connections between the approaches 

Since from (^Sf) the structure functions are real, one can rewrite (^) as 

oo 

m=—oo 

oo 

m=— oo 

Replacing mhj m + k in the first term of this sum and using (p^Bf) then gives 



an alternate formula for the average probability distribution in ( ^3]) 



P{xo + kVt,t) = (4L)-^ 



Y k™ - ^^+fce^(2'"+*^)-^o/^ -2|V^fc/2 



(24) 

The above expression clearly indicates a link between the average proba- 
bility distribution and interference of the amplitudes ipm and ipm+k- Further, 
this interference is seen to be destructive when xq = IL for some integer / 
with kl even, thus recovering condition (||) for channel structures. The de- 
structive interference approach is thus a special case of the travelling wave 
approach, where condition (|1]) corresponds to particular minima of the struc- 
ture functions (at z = I) for the case of slowly varying energy amplitudes. 
Finally, ([T^) can be used to rewrite (03) in terms of energy amplitudes 
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ipm with m > 0: 



\k\-i 



P{kVt + IL, t) = L-^Re < 1 - E ^*m^m+\k\ + ^ E ^*m^\k\->n 

[ m=\ ^ m=l (m^|fc|/2) 

(25) 
when hi is an even integer. Thus the approximation in (^) is seen to be quite 
good for small \k\^ and indeed is exact for \k\ equal to 1 or 2. 



4 Properties of structure functions 

4.1 Wigner function form 

Expression (^) for the structure function S^ is not always convenient to 
use. It generally involves an infinite summation, and the contributing energy 
amplitudes are not always easily calculated. It would therefore be useful to 
have a formula for 5*^ which is directly related to the initial wavefunction 
V'(x,0). 

To obtain such a formula, first define the normalised wavefunction 0(a;) 
on the interval (— i^, V) by 

0(x) = 2-^l\^{x,^)-^{-x,^)\ (26) 

oo 

= L-i/V(2z) E ^ne^'^-^•/^ (27) 



where the second line follows via (p!7D and (pTSl) . Thus </)(x) is an antisym- 
metric extension of the initial wavefunction ^(x, 0). Now let iy<^(a;,p) denote 
the Wigner function of 0(x) [[lOl : 

/oo 
t/#*(x-l/)0(x + l/)e-2*P^/^ (28) 

-oo 

For < a; < L one then has, as shown further below, the remarkable relation 

5,(x/L) = -Kh [iy^(a;,pfc) + {-VfW^ix - L,p,)] - l^fc/2|V2, (29) 

connecting 5*^ with W^, where 

Pfc = Tink/{2L) = kMV. (30) 

13 



This expression can be trivially extended to evaluate Sk{z) for all values of 



z via the first of the relations in (22). 

As will be seen below, ( P^j ) provides a very convenient method for evalu- 
ating structure functions, and also for analysing the dependence of quantum 
carpets on various properties of the initial wavefunction. Expansions of the 
probability distribution P{x, t) in terms of Wigner functions have been pre- 
viously obtained @, 11]. However, these expansions lead only to expressions 



for Sk involving infinite sums of Wigner functions, in contrast to (pQ]). The 
latter may be derived by substituting ( p6|) and (|30|) into ( p8|) and recalling 
that 0(x) vanishes for |x| > L hj definition, thus yielding 



ATrhLW^{x,pk) 



/Mix) 
^y^i{m+n-k).y/L 

+ 2 Y. C^fc-me-*('™-')""/^M(a;), 



where 

M[x) := max{L — |x|, 0}. 

Substitution into the righthand side of (|29|) gives (0) as required. 



(31) 



4.2 Localisation and symmetry effects 

To see how localisation and symmetry properties of the initial wavefunction 



directly affect the structure functions, one may substitute (|26|) and (pSD into 
to obtain 



Su{xlV) = {nh/2)[w4x,pk) + i-l)''W4L-x,-pk) 
~{l/2)\l4x,pk) + \^k/2\^' 



(32) 



for < a; < L, where W^ denotes the Wigner function of ip{x, 0), and the 
"interference" term I^{x,pk) is given by 



mx.Pk) 



IZo dyi'*{y + 2;, 0)^(2/ - x, o)e'^''y/^ + c.c, o<x<l/2 



-iriiL-x,-pk) 



L/2<x< L. 

(33) 



14 



Now, for example, suppose that the initial wavefunction ip{x, 0) is well- 
localised about some point Xq 7^ 1/2. The corresponding Wigner function 
W^{x,p) will then be similarly localised, and hence from (^) and ( P^D one 
predicts (i) a ridge structure associated with the trajectory xq + kVt; and 
(ii) a ridge/channel structure associated with the trajectory L — xq + kVt 
for even/odd values of k. An example verifying this prediction is given in 
section 5.2 below. 

It is also of interest to consider the cases of trajectories of the form x = 
IL + kVt (arising from the destructive interference approach in section 2). 
From m), (m) and flSSD one has 



^,(0 = (-i)"'5,(o) 



' _'^\kl+l 



dxP{x, 0) cos{TTkx/L) + \ipk/2\ /2 




(34) 



where P{x, 0) is the initial position probability distribution. The average 
probability density along these trajectories thus only depends on the phase 
structure of iplx, 0) via ipk/2 (which is zero for odd values of k and typically 
small in general). 

If initial distribution P{x, 0) is well-localised within a region [xq — Sx, Xo + 
6x) such that 

6x/L<^\k\-\ (35) 

then from ( [5^ ) the average probability density along x = IL + kVt is essen- 
tially determined by the properties of the cosine function cos nkxa/L. For 
example, if Xq satisfies 

kxo/L = j (36) 

for some integer j, then from ( ^5] ) and ( |5^ ) one has 



P{lL + kVt,t)^{l-{-iy+''')/L, (37) 

corresponding to channel and ridge structures as j + kl is odd and even re- 
spectively (the ^fc/2 term has been ignored in (|37D, as from ([TT|) it is typically 
negligible for locahsed wavefunctions). Similarly, if j is replaced by j + 1/2 
in ( |36D then the (—1)'^^+-' term in (|37| ) vanishes, and any linear structure 
associated with the trajectory is suppressed. 
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Finally, rather than supposing the initial distribution to be well-localised, 
consider instead the case where P{x, 0) has an approximate reflection sym- 
metry about some point x*, i.e., 

P{x* + x,0)^P{x* -x,0). (38) 

Noting that P{x, 0) vanishes outside (0, L) one then has from (p^ that 

/•oo 

{-l)''^Sk{l) ^ - dxP{x* + x,0)\cos7ik{x* + x)/L 
Jo 

+ cos7rA;(x* — x)/L] — \^k/2\'^ /2 

= -2cos^— / rfxP X* + X, cos ^ - -hAfc/2^. 39 
L Jo L 2 

Thus the structure function is modulated by the cosine function cos7ikx*/L. 
Linear structures associated with the trajectories x = kVt + IL are therefore 
enhanced when kx*/L = j for some integer j, and suppressed when kx*/L = 
j+1/2. Examples of such enhancement /suppression are given in the following 
section. 

5 Examples 

5.1 Uniform initial wavefunction 

The case ip{x, 0) = L^^/^ on (0,L) was considered by Berry 0, who showed 
that the corresponding probability landscape was a fractal and explained 
the observed channel structures via destructive interference (Example (ii) of 
section 2.3 above). For this case it follows from (0) and (|18D that 

iJn = 2^2/ (nvr) (40) 

for odd n, with ■?/'„ vanishing for even n. The structure function S2k{z) can 
be evaluated via any of (^TD, (p9|) or (^2]), to give 



S2k{z) = 2{nk)-' sin[27rfcmin{2, 1 - z}] - 4o - \i^k\V2 (41) 

for < 2; < 1, where 6m is the Kronecker delta and sin[x]/x is evaluated as 
1 for a; = 0. As usual this may be extended to other values of z via 
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In Figure 4 (^) is plotted for k = 1,2 and 3. It is seen that while channels 
are associated with trajectories x = 2kVt + lL/2 for even values of (/c + 1)/, as 
predicted in Example (ii) of section 2.3, these channels are relatively broad 
and shallow. Indeed the average probability density along these channels 
follow from (^2]), (^) and (^) as just the background probability 1/L for 
even k, and [1 — 4/(7rA;)^]/L for odd k. Channel visibility thus decreases 
rapidly as |fc| increases. 

Since the initial probability distribution is symmetric about x* = L/2 for 
this example, one expects from ( ^9|) that channels with slopes of {2k + 1)V 
are suppressed. Indeed, as is most easily seen from ( p^ (recalling ipn vanishes 
for even values of n), one finds 5*2^+1 (^) = 0. Hence no linear structures of 
slope {2k + 1)V are predicted, and ( pil| ) and (^TJ) then yield the surprisingly 
simple decomposition 

CO 

P{x,t) = 2 Yl {nkL)-^(r{x-kVt)sm[27rk{x-kVt)/L] + [l-2a{x)]/L (42) 

k = — OD 

of the fractal probability landscape, where a{x) is defined to be +1, 0, and 
—1 as the fractional part of x/L is respectively less than, equal to, or greater 
than 1/2. 

5.2 Gaussian initial wavefunctions 

Consider now the inital wavefunction 

^(x, 0) = ir(27ra2)-i/^e-("-^)'/(^'^') e'^", (43) 

where < x < L and i^ is a normalisation constant. It will be assumed that 
0" -C a;, L — x, i.e., that the wavefunction is well-localised on the interval. One 
may then make the extremely good approximation i^ = 1, so that ip{x, 0) is 
effectively a Gaussian wavepacket centred at x with average momentum p. 
It is easiest to evaluate the structure functions for this case via (B2D and 



(|33D , where one makes the (again extremely good) approximation that (^) 
can be extended over the entire x-axis. Performing the resulting Gaussian 
integrals then yields 



\LJ 2 



-2<T^(pfc-p)^ -(L-x-x)'^ -2<T^(pfc+p)^ 



;-i)'[e- 
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cos^{pkX-px), 0<a;<| 

2i2 COS ^{pkX + p{L — X)), -^ < X < L. 

The term — (l/2)|'?/'fc/2p in ( ^21) has been ignored, as it vanishes for odd k 
and from ([TD|) and (|I2]) is only of order cr/L ^ 1 for even k. 

The structure function 5*1(2;) is plotted in Figure 5 with x = L/3 and 
(7 = L/40, for the cases p = (dashed line) and p = Ihnfi/L (solid line). 
In the latter case it is seen that there is a deep channel associated with 
the trajectory x = Vt, flanked by a high ridge on the left and a moderate 
ridge on the right. Further, from (p^) one has 5'_i(l + z) = Si{l — z), 
and hence there is, conversely, a high ridge associated with the trajectory 
X = L — Vt, flanked by a deep channel on the right and a moderate channel 
on the left. These predicted features following from Figure 5 may be directly 
observed in Figure 6, where the probability landscape corresponding to initial 
wavefunction (|43|) is plotted for x = L/3 and p = Ibirh/L. The structure 
functions thus accurately predict the shapes of the observed linear structures. 

For the p = case (dashed line) in Figure 5 it is seen that as well as 
channel and ridge structures at the endpoints, there is a ridge associated 
with the trajectory x = L/3 + Vt, and a channel associated with the tra- 
jectory X = 2L/3 + Vt. These latter structures correspond to the Wigner 
functions W^{xo,pk) and W^{L — XQ,pk) discussed in section 4.2 for localised 
wavef unctions, and can again be observed in a plot of the probability land- 
scape. The corresponding structures for the p = Ibnh/L case are observed to 
be much less pronounced. This is because W^{x,p) is localised about p = p, 
and is thus relatively small at p = pi. 

Gaussian initial wavefunctions also provide an example of the prediction 
in ( pTl) for well localised initial wavefunctions. In particular, from (^^ (re- 
calling the assumption a <t^x, L — x), one finds that 

<S'fc(0) '^ -e 2l2 cosnkx/L. (45) 

Hence (^) follows whenever conditions ( P^D and ( |5BD are satisfied (with 
6x = a and xq = x). Further, channels and ridges associated with trajectories 
X = kVt + IL are seen to be suppressed if nkx/ L = j + 1/2 for some integer 
j. Since Gaussian initial wavefunctions are not only localised about x, but 
are approximately symmetric about x, one may also obtain (|45|) from (|39|) , 
with X* = X. 



6 Periodic gratings and rigid rotators 

Consider now a plane wave of wavelength A incident on a one- dimensional pe- 
riodic grating, of period 2L in the x direction. In the paraxial approximation 
the amplitude of the diffracted light then has the general form 



{x,z} 



(2L)-l/2 ^ ^^^i{nx-n^Vz)./L^ (45) 



where z measures distance propagated perpendicularly to the grating, the 
Fourier coefficients 0„ are determined by the initial diffracted wave 0(x, 0), 
and1/ = A/(4L). 



Clearly p6|) is formally similar to (|T^ for a quantum particle moving 
freely between two endwalls, and indeed analogues of quantum carpets in 
the grating context have been previously observed and investigated |0, |]. 
The significant differences between (^) and (|l^) are (i) 0(x, z) is periodic 
on the entire x axis whereas ilj{x,t) vanishes outside the interval (0,L); and 
(ii) the coefficients 0„ need not be antisymmetric as per (plsl) . 

Many of the results of sections 3 and 4 may be directly translated into 
the grating context. For example, the light intensity distribution I{x, z) = 
\(f){x,z)\'^ has the travelling wave decomposition 



I{x,z) = {2Ly 



k=—oo \ / 



(47) 



in analogy to (OT), where Iq is the integrated light intensity per period of 
the grating and Sk is given by (|2TD with ^„ replaced by v^0„. Similarly, the 
analogue of (ESI) is 



J(xo + kVz.z) = (2L)-i[/o + Sk{x^/L)l (48) 

and thus the structure function S^ determines the average locations, shapes, 
etc of relatively dark and bright structures of slope kV in the x-z plane. The 
average intensities in (^) are precisely the "lane contrast functions" first 
defined in [|] for predicting the locations of dark structures. 

While one may also write down a result analogous to (^Uj) , an alternative 
Wigner function expression for the structure functions is more useful in the 
grating context. In particular, define the periodic Wigner function 

W^{x,v)=T^-' r rf#*(x-y,O)0(x + y,O)e-2*P^. (49) 
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This differs from ( |28| ) in that does not vanish outside the interval {—L, L). 



One may then show, similarly to the proof of (29), that 

Sk{x/L) = nW^{x, 7Tk/{2L)) - |4/2p, (50) 

where (f)k/2 is defined to vanish for odd k. 

As an example of (0), suppose that a plane wave is incident on a sinu- 
soidal phase grating, as studied in |^, i.e., 

^{X, 0) = (/o/2L)l/2 ^^acos{nx/L)_ (5^^ 

Substitution into (^) and using the standard Bessel integral 

/ de cos[ke - a sin 9] = vr Jfc(a) (52) 

Jo 



then yields, via (EHf) and (|50|) 



/(xo + kVz,z) = /o/(2L) [1 + {-lfJki2asm[nxo/L]) - |</)fc/2rj . (53) 

Here \(f)k/2\ is zero for odd k, and equal to Jk/2{(^) for even k |Q]. 

The last term of (0) provides a correction to equation (25) of [||, where 
the average intensity along the trajectory x = Xo + kVz was calculated for the 
same example. More generally, (|50|) provides a convenient way to calculate 
the structure function directly from the initial amplitude (f){x,0), and (^) 
provides a simple decomposition of the intensity distribution as a sum of 
travelling waves. 

The above results may be immediately translated into the context of the 
two-dimensional quantum rigid rotator, with Hamiltonian H = J-^ /{2I). In 
particular, 0(^, t) may be interpreted as the time- dependent phase ampli- 
tude of the rotator, where one takes L = n and V = h /(2/). The carpet 
structures for this case lie on the cylinder generated by the phase and time 
coordinates 6 and t. 



7 Conclusions 

The carpet structure of a quantum particle confined between two endwalls 
has been greatly elucidated by the travelling wave decomposition (^) of the 
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probability density, and the details afforded by this decomposition are seen 
to go well beyond the ambit of the destructive interference approach. In par- 
ticular, the structure functions Sk in (^^ provide all information pertaining 
to the average properties of linear structures in the probability landscape. 

The structure functions may be evaluated either from the energy am- 
plitudes of the wavefunction, via (pTI); or from the initial wavefunction, via 



(|32D . The latter formula permits the dependence of generic features of carpet 
patterns to be directly determined from corresponding properties of V^(x, 0). 
Similarly, for the periodic grating and quantum rigid rotator one may use 
( |50D to directly evaluate structure functions from 0(x, 0). 

One advantage remaining to the destructive interference approach is that 
it may be generalised to predict the carpet structure of semiclassical wavepack- 
ets moving in arbitrary one-dimensional potentials (see Appendix). It is 
hoped that a corresponding generalisation of the travelling wave decomposi- 
tion ( pOD can be found. 
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Appendix 

It will be briefly indicated here how the destructive interference approach 
of section 2 may be generalised to predict the curved carpet structures arising 
in the probability landscapes of one-dimensional particles moving in general 
potentials m^, |] . For semi-classical wavepackets the analysis is very similar 
to that of subsection 2.1, and the results agree with those of Kaplan et al 
obtained by consideration of degeneracies of the propagator for P{x,t). 

Now, for high energies the n-th energy eigenstate of a particle moving 
in a one- dimensional potential U{x) with two classical turning points is well 
approximated as |T2| 



il^nix) = 2[M/Tr,Pn{x)Y/^sm[h-' f Pn{x)dx + 7r/4], (54) 

where M is the particle mass; T„, p„(x) and Xn are the period, momentum 
and leftmost turning point respectively of a corresponding classical orbit of 
energy En'-, and En is implicitly deflned by the Bohr-Sommerfeld rule 

Pn{x)dx = {n + l/2)2Tih. (55) 

The appearance of the sine function in (|5^ is analogous to that in (|^), 
and for an initial superposition ip{x,Q) = I]c„'?/'„(x) of such states one can 
easily obtain the analog of @: 

oo 

^P{x,t) = 5:[A„,+,fc|(x)e^<^±(^'*'"+l^l)-A„(x)e^<^^(^'*'")" 

n=l 

\k\ 

+ 5]A„(x)e^'^±(^'*'"), (56) 

n=l 

where 

An{x) = [M/TnPn{x)Y'^Cn. (57) 

0±(x,t,n) = ±h'^ f pnix)dx-n'^Ej±TT/A. (58) 

As in (H), the upper (lower) phase subscript is chosen in (|56D when k is 
positive (negative), and an overall phase factor has been dropped. 

Destructive interference in the flrst summation in (^), giving rise to a 
channel structure, can take place along a spacetime trajectory Xk{t) if the 
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phase condition (^ holds for (l)±{x,t,n) defined in (|58| ) (providing that the 
amphtude functions An{x) in (|57D vary sufficiently slowly with n). Substi- 
tuting ( pBD into (|^) and differentiating with respect to t yields the condition 

dxkjt) En+\k\ — En , . 

dt Pn+\k\{Xk) + Pn{Xk) 

for the trajectory Xk{t), where the +/— sign is chosen according to whether 
k is positive/negative. Moreover, if the trajectory passes though point Xq at 
time to, one requires from (^ and (|58|) that 



(En+k - En)to - / Pn{x)dx - / Pn+k{x)dx 



/h = 7r/2 mod 27r. (60) 



For the case of a potential energy which is symmetric about xq, this expres- 
sion may be simplified via (|55| ) to give 



{En+k - En)to/h = (n + k/2 + 1)71 mod 27r. (61) 

Channel structures for various potentials have been numerically observed 
for a number of examples [^ |^, Q , and (|59| ) has been previously derived in 
[§], via a decomposition of the probability distribution rather than of ^(x, t) 
(which produces an extra set of trajectories, discarded as "classical"). Here 
it is seen that (^) arises directly from destructive interference of energy 
amplitudes, in a manner entirely analogous to the case of the particle in a 
one-dimensional box. 

It is hoped to further investigate conditions ( ^9|) and (|60|) elsewhere. 
Here a simple prediction generated by (|59D will be pointed out. In par- 
ticular, if the energy eigenvalues En increase slowly over ranges of length 
\k\ (at least for values of n for which the amplitudes An{x) are signifi- 
cant), and Pn{x) varies slowly over such ranges, then from (|59|) one has 
dxk/dt ~ k{dEn/ dn) / {2pn{xk)) . Thus, if two channels coresponding to two 
values k and k' intersect at some point in spacetime, their slopes at the point 
of intersection are predicted to be approximately in the ratio 

{dxk/dt)/{dxk'/dt) ^ k/k'. (62) 

This provides a simple test of the applicability of this approach to a given 
quantum carpet structure: the strongest channels, corresponding to small 
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values of \k\, are predicted to intersect with slopes related by simple ratio- 
nal numbers. Note from (|56|) that conditions (^) and (|60| ) for destructive 
interference need in fact only hold over the range of n for which the ampli- 
tudes An{x) are significant, to ensure effective destructive interference. This 
is fortunate, as these conditions cannot in general hold for all n; however, it 
implies that this approach can in general only be applicable to superpositions 
of a relatively narrow band of energy eigenstates. 
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FIGURE CAPTIONS 

Figure 1. Density plot of (part of) tlie probability landscape for the 
initial wavefunction in (JT^, with A^ = 20. The dark channel observed to 
run from the bottom righthand corner to the top lefthand corner of the plot 
corresponds to the trajectory x = L/3 — 3Vt, and is not predicted by the 
destructive interference approach. 

Figure 2. Plot of 1+Ss{z) for the initial wavefunction in (|16|) , with A^ = 
20. From (^), the sharp minimum in the vicinity oi z = 1/3 corresponds to 
a channel in the probability landscape, along the trajectory x = L/3 — 3Vt 
(as observed in Figure 1). 

Figure 3. Approximate reconstruction of the density plot in Figure 1 
via the travelling wave decomposition (PDj), where only terms with |/i;| < 5 
have been included. All linear structures of slope less than or equal to 5V in 
magnitude are successfully reproduced. 

Figure 4. The structure functions 82(2) (solid line), S4{z) (dot-dashed 
line) and 5*6 (-2) (dotted line) for the case of a uniform initial wavefunction, 
plotted via (^). 

Figure 5. The structure function Si{z) for the approximate Gaussian 
initial wavefunction in (^3]) with x = L/3 and a = L/AO, for the cases p = 
(dashed line) and p = Ibrch/L (solid line). 

Figure 6. Density plot of the probability landscape for the approxi- 
mate Gaussian initial wavefunction in (|43| ) with x = L/3, a = L/40 and 
p = 15iTh/L. The detailed shape of the linear structure running from the 
bottom lefthand corner to the top righthand corner corresponds to that of 
the structure function of Figure 5 (solid line) in the neighbourhood of 2; = 0. 



26 



L/9V 




L/3 



2.0 



1.5 



1.0 



0.5 




L/9V 




L/3 



0.3 . 



0.2 



0.1 



0.0 



-0.1 



-0.2 



-0.3 . 



-0.4 




1.0 [ 



0.5 



0.0 




-0.5 



\ 

\ 



/ \ 

/ \ 

/ \ 

/ \ 



-1.0 







1/3 



2/3 



1 



L/V 




